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Abstract 

We use the AdS / CFT correspondence to perform a numerical study of a phase transition 
in strongly-coupled large-A^^^ A/" = 4 Super- Yang-Mills theory on a 3-sphere coupled to a 
finite number Nf of massive J\f = 2 hypermultiplets in the fundamental representation of 
the gauge group. The gravity dual system is a number Nf of probe D7-branes embedded in 
AdS^ X S^. We draw the phase diagram for this theory in the plane of hypermultiplet mass 
versus temperature and identify for temperatures above the Hawking-Page deconfinement 
temperature a first-order phase transition line across which the chiral condensate jumps 
discontinuously 
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1 Introduction 



The AdS/CFT correspondence [1] equates the low-energy effective theory of string theory, 
supergravity (SUGRA), on the background AdS^ x with strongly-coupled = 4 SU{Nc) 
super- Yang-Mills theory (SYM) in the large- A^c hmit hving, in some sense, on the bound- 
ary of AdS^. At finite temperature this four-dimensional boundary is x S^. This field 
theory contains only adjoint fields and exhibits a first-order deconfinement transition at the 
Hawking-Page temperature, corresponding to AdS-Schwarzschild black hole condensation 
[2]. The thermodynamics of this SYM theory has been studied in great detail both at strong 
coupling via AdS / CFT and at zero [3] and weak coupling [4] . 

Why consider a theory on a compact space? For this conformal theory, the deconfinement 
transition cannot occur in infinite volume, i.e. on spacetime x because with no scale 
nothing can set a transition temperature. The introduces a new scale, and indeed the 
deconfinement transition occurs at a temperature of order the inverse radius [4]. On 
the gravity side, this is the statement that the Hawking-Pagc transition is only apparent 
in global AdS coordinates. Poincare patch coordinates are limited to the high-temperature 
black hole solution, corresponding to the deconfined phase in the SYM theory. Of course, no 
phase transitions can occur with a finite number of degrees of freedom, which then motivates 
the 't Hooft limit. 

In large- A'c Yang-Mills, or more generally in gauge theories with only adjoint fields, the 
deconfinement transition comes from order dynamics. One way to identify the decon- 
finement transition is to ask whether the order A^^ contribution to the free energy is zero 
or not. In pure Yang-Mills, for example, at sufficiently low temperature this contribution to 
the free energy will be zero, as the degrees of freedom are color singlet glueballs of whom 
there are order (in the large- A'c book-keeping), while at sufficiently high temperature it 
will be nonzero, indicating that the contributing degrees are freedom are deconfined gluons 
of which there are order N^. This is what the word "deconfinement" means in the M — ^ 
SYM theory on a 3-sphere studied at strong coupling using AdS/CFT. 

Introducing fundamental fields into this field theory requires introducing an open string 
sector. In reference [5] D7 flavor branes were introduced to the usual AdS/CFT construction, 
which equates the field theory on a stack of coincident D3-branes with gravity in the geometry 
supported by the 3-branes. Keeping the number Nf of D7's finite while taking the number 
Nc of D3's to infinity will produce in the near-horizon limit probe D7-brancs embedded in 
AdS^ X S^. These branes will wrap the AdS^ and an 5'^ inside the factor. The probe 
branes break half the supersymmetry. 

In the field theory this corresponds to adding to the J\f — A theory a number Nf olM — 2 
hypermultiplet fields (henceforth "hypers") in the fundamental representation of SU{Nc). 
The theory will remain conformal in this Nj/Nf. limit. The hypers can be given a mass 
by "ending" the D 7- branes at some finite value of the AdS radius [5]. Such masses obviously 
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break the conformality but will preserve the M = 2 supcrsymmetry. On the gravity side this 
is manifest as super symmetric probe D7's ending at finite radial coordinate, which clearly 
breaks the radial isometry. 

The addition of fundamental fields introduces the possibility of a chiral phase transition 
analogous to that in quantum chromodynamics (QCD) or at least large-iVc QCD. By the 
usual large- A^c counting rules, all dynamics of fundamental fields is suppressed by a factor 
of l/Nc- For large- A^c theories, then, the chiral transition involving fundamental fields is 
in some sense a small perturbation, suppressed by a power of on top of the order N"^ 
dynamics of the deconfinement transition. 

Most of the remainder of this section, as well as Sections 2 and 4, contains little new 
material and is presented to put our work in the proper context. Readers familiar with the 
topic may want to go directly to Sections 3 and 5, where we present our methods and new 
results. Our main result is figure (2), as summarized in the last three paragraphs of this 
introduction. 

We now review the chiral phase transition in QCD, to compare and contrast with our 
H — 2 supersymmetric theory on a compact space. We will be careful about two limits: 
large- A^c and zero quark mass. In general (any and A^/), for massless quarks chiral 
U{Nf)L X U{Nf)R ^ SU{Nf)L X SU{Nf)R x U{1)b x U{1)a is an exact symmetry of the 
Lagrangian. At zero temperature, SU{Nf)L x SU{Nf)fi is spontaneously broken to the 
diagonal SU{Nf)v- The chiral symmetry is restored at sufficiently high temperatures. The 
order parameter for the chiral transition is the quark bilinear condensate, which is nonzero at 
low temperatures and zero at high temperatures, signaling the chiral symmetry restoration. 
Quark masses explicitly break the chiral symmetry, but we will follow the convention of 
calling the quark bilinear condensate an order parameter nonetheless. 

In QCD with massless quarks, the order of the chiral transition depends on Nc and Nf. 
We start with A"c = 3. For Nf > 3 massless quarks, universality arguments indicate that 
the chiral transition is first order [6]. For Nf = 2, the chiral transition is believed to be 
second order [7, 8]^. For Nf = 1, things are more subtle. Now the symmetry is simply 
U{1)b X U{1)a and the question becomes whether the axial symmetry, which is anomalous 
at T = 0, is restored as the temperature rises. The instantons responsible for the anomaly 
are suppressed as T rises, so in fact the axial symmetry is believed to be restored at least 
in the strict T = oo limit. As the decrease in the density of instantons is believed to be 
smooth, no chiral transition exists for Nf = 1 [6] as a function of T. Turning to the large- A'c 
limit, for Nf > 2 the theory has no stable IR fixed point and hence the transition must be 
first order [6]. At large A^c; the axial anomaly is suppressed and to our knowledge it is not 
known whether there is a phase transition for Nf — 1 and if so at what order. For the theory 

^Lattice simulations have not settled the issue of the order of the transition for Nf = 2. As just a small 
sample: [9] found results consistent with a second order transition while more recently [10] found indications 
of a first order transition 
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we are studying the relevant chiral symmetry is also U{1)b x U{1)a, so it is precisely this 
Nf = 1, large Nc case that is closest to the supersymmetric theory we have control over. The 
same symmetry also governs a single staggered fermion in the strong coupling limit. Lattice 
results for that theory [11] indicate a second order transition. 

Next we wish give the quarks nonzero masses and ask what happens to the deconfinement 
and chiral phase transitions as we dial the mass. Dialing quark masses is not possible in 
nature but is possible on the lattice, hence we turn to lattice results although we will only 
present them schematically. The phase diagram for Nc = 3 and Nf = 2 + 1 , is shown in fig. 
1 (a.) [12]. Here 2+1 means the quark masses are held at fixed ratios m = m„ = ma = riris 
for some r < 1. What appears are two lines of first-order phase transitions, each ending in 
a critical point. The lower line is that of chiral transitions, ending in the point C; the upper 
is the line of deconfinement transitions, ending in the point D. 

The chiral transition is identified by a discontinuous jump in the chiral condensate. Only 
on the m = axis, where the chiral symmetry really is a symmetry, is this a genuine 
symmetry breaking transition. 

The deconfinement transition is identified by a discontinuous jump in the Polyakov loop, 
which is the time-ordered exponential of the holonomy of the gauge field around the time 
circle. This measures the response of the system to an infinitely massive source of color 
charge (a quark), and may be written 



where AF is the difference in free energy between two equilibrium states, one with an 
infinitely massive test quark and one without. For pure glue in the confined phase this 
difference is infinite and hence {Vt) — 0. One can think of this as the statement that the 
color flux lines from the test quark have no place to end except inflnity and hence because the 
flux lines have finite energy density the total energy is infinite. In the deconfined phase, the 
fiux lines are screened, hence (Vt) is finite. Adding dynamical quarks, massive or otherwise, 
changes this since now flux lines from the test quark have someplace to end. In other words, 
with dynamical quarks (Q) is always nonzero, although it may still jump discontinuously. 
Only in the m — > oo limit where the quarks decouple is the pure glue story recovered. 

In summary, both order parameters are nonzero everywhere on the interior of fig. 1 (a.). 
At m = oo where the quarks decouple, the deconfinement transition occurs at roughly 
Tfi ~ A.QCD ~ 175MeV. Similarly, for m = where the chiral symmetry is an exact 
symmetry of the Lagrangian, the chiral transition must occur at the same scale (that then 
being the only scale in the theory) . The horizontal line of realistic quark masses falls between 
points C and D, hence a smooth crossover rather than a phase transition is expected in the 
real world. 

Now we consider the same diagram but with N^ oo, Nf finite. The answer, again 
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Figure 1: (a.) Phase Diagram of QCD with N^^ 3, Nf ^2 + 1 (b.) N^^ oc, Nf ^2 + 1 

only quahtative, is depicted in fig. 1 (b.). The deconfinement transition, being order N^, 
does not care about the quark masses and is simply a vertical line. We have drawn the 
chiral transition line the same way although in principle this need not be the case. Large- A^c 
arguments indicate that the chiral condensate is independent of T in the confining phase 
[13, 14] and hence the chiral transition must occur at a temperature T^h > 7d, we have thus 
drawn the deconfinement line to the left of the chiral transition. 

In this paper we study the chiral transition in a theory that is in some sense simpler than 
QCD: large- A^c M — 2 SYM on a 3-sphere. At zero temperature and zero mass, conformal 
invariance precludes a nonzero chiral condensate, but finite temperature and finite mass both 
break conformal invariance. We will draw the phase diagram for this theory in the plane of 
hyper mass versus temperature, identifying the line of a first-order chiral phase transition, 
which occurs only in the deconfined phase, i.e. above the Hawking-Page temperature. This 
transition is, in strict thermodynamic terms, not a chiral symmetry breaking transition, 
but the chiral condensate does jump by a finite amount so for simplicity we will refer to 
it as a chiral transition. In contrast to QCD we will see that chiral symmetry is unbroken 
for zero mass at all temperatures, so the deconfinement transition does not come with a 
standard chiral symmetry breaking transition. Another key difference from QCD will be an 
independence on the number Nf. Where the order of the chiral transition in QCD depended 
on Nf, in our theory it is first order for any finite Nf such that Nj/N^ 0. The J\f = 2 
Lagrangian explicitly breaks the non-abelian chiral symmetry down to its vector subgroup 
due to a Yukawa coupling, so the chiral symmetry breaking we are analyzing is a breaking 



of U{1)a X U{1)b down to U{1)b as in iVj = 1 QCD. 
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The thermodynamics of large- A^c M — 4, SYM on a 3-sphcrc at zero 't Hooft couphng, 
but with a Gauss' law constraint requiring color singlet states, was studied by Sundborg [3], 
who found that the deconfinement transition was first order. Aharony et. al. [4] studied the 
same theory at small but finite coupling and found that the deconfinement transition could 
be either a single first-order transition or two continuous transitions. For pure Yang- Mills on 
a 3-sphere, the former scenario was shown to be the right one [15]. Schnitzer [16] introduced 
flavor but with Nf/Nc flnite, flnding that the deconflnement transition becomes third order. 

The chiral transition of the D3-D7 system for the Af = 2 theory in flat space was studied 
in [17, 18]. A first-order transition was found [18, 19]^. A similar phase transition was also 
observed for fiavor branes in the supergravity dual of confining gauge theories [20]. The 
phase diagram must by conformality be linear: the quark mass is the only scale that can 
set a transition temperature. As we will argue below, this infinite volume result can be 
interpreted as a high-temperature result for the theory on the 3-sphere, and indeed we will 
find a straight line at high temperatures, providing a useful check of our methods. 

Our final result for the chiral transition is the phase diagram in figure (2). The vertical 
line is the Hawking-Page deconfinement temperature, Tfjp = 3/27r (in units where the AdS 
radius is one), which comes from order A*"^ dynamics and is independent of the quark mass. 
Below Thp, the fiavor brane action is independent of temperature, hence whatever occurs 
at T = can be extended up to Tffp. This is consistent with the large- A^c field theory 
arguments that the chiral condensate is independent of temperature in the confining phase 
[13, 14]. We find no first-order phase transition as a function of temperature below Thp- 
This comes as a surprise since at zero temperature we see a topology change for the D7 
brane as a function of mass. 

Above Thp, the flavor brane action depends on the temperature via the presence of the 
black hole horizon. In this temperature regime, the chiral condensate and free energy exhibit 
behavior characteristic of a first order transition as m increases, i.e. the first derivative of the 
free energy is discontinuous and so on. We have drawn the first-order chiral transition line. 
Below Thp, the chiral condensate smoothly changes, as does the free energy, as m — > oo, so 
we draw no first-order transition line. Everywhere on the m — axis, we find numerically 
that the chiral condensate is zero. 

This paper is organized as follows. In section 2 we review the inclusion of flavor branes into 
the Af — 4 theory and identify the chiral symmetry in this theory and its order parameter. 

^Whilc wc were finishing tliis manuscript two papers appeared which also reanalyzed the phase transition 
in this theory, [21] and [22]. In both papers the picture of a first order phase transition as found in [18, 19] 
is confirnied. Like us, they simplify the numerics by integrating outward instead of shooting inward. In 
addition [21] argues that when the flavor brane touches the horizon the probe action acquires a new scaling 
symmetry and a self-similar structure. Their general analysis of Dp/Dq systems corresponding to large- A^c) 
(p-l-l)-dimensional thermal gauge theory in infinite volume with fundamental matter confined to a (q-|-l)- 
dimensional defect also shows that such phase transitions must be first order. Our results are an extension 
of this since we consider finite volume and find again a first-order transition in the deconfined phase. 
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Figure 2: The phase diagram for flavored Af — A Super- Yang Mills on a sphere as a function 
of temperature and hyper mass. 

In section 3 we write down the probe branc action and counterterms, explain how we extract 
the condensate and free energy using holographic renormalization, and explain the boundary 
conditions for our numerics. In section 4 we reproduce known results in infinite volume. In 
section 5 we present numerical results for the embeddings, condensate, free energy and finally 
for the phase diagram in the finite volume case. We conclude in section 6 with possibilities 
for extensions of this study. 

2 Adding Flavor to AdS/CFT 

We will now describe our system in greater detail, clarifying both the D-branc construction 
and the symmetries of the boundary theory. The AdS / CFT construction starts with a stack 
of coincident D3-branes in fiat ten-dimensional space in the limit where the number of D3- 
branes, N^, goes to infinity and the spacetime in the near-horizon limit becomes AdS^ x S^. 
The limit of zero string length a' — > then decouples open and closed strings, resulting in the 
now-standard correspondence between the low-energy effective theory of closed strings, ten- 
dimensional supergravity (SUGRA) in this background, and the low-energy effective theory 
of open string modes on the D3 worldvolume, Af = A super- Yang-Mills (SYM) theory in four 
dimensions in the large- Ai'c, strong 't Hooft coupling regime. This four- dimensional conformal 
field theory (CFT) will "live" in a spacetime with the topology of the AdS;^ boundary. Our 
interest here is in global thermal AdS (and Ad^'-Schwarzschild) coordinates, with boundary 
topology X S^, which we will call the "curved case", in distinction to the infinite- volume 
"flat case", that is, Poincare patch coordinates with boundary topology x R^. We will 
use the limit of infinite volume to check our methods against known fiat case answers. 

Thermal AdS is beheved to undergo a first-order phase transition [23, 2] called the 
Hawking-Page transition. At low temperatures the picture is of thermal radiation in equilib- 
rium with itself (assuming perfectly reflecting boundary conditions at the AdS boundary), 
whereas when the temperature rises the spacetime undergoes a topology change to the AdS- 
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Schwarzschild solution, or in other words a black hole condenses. This transition is only 
apparent in global AdS coordinates, while the flat case is strictly in the high-temperature 
black hole phase. 

In the dual J\f = 4 SYM boundary theory, the Hawking-Pagc transition corresponds to a 
large-iVc deconfinement transition. "Deconfinement" here means the spontaneous breaking 
of the center symmetry and the associated behavior of that symmetry's order parameter, 
the Polyakov loop, which jumps discontinuously (from zero to nonzero) as the temperature 
rises. This transition is only apparent for the theory on a compact space, where the inverse 
radius can set a transition temperature. 

In fact, the ratio of the and radii is the only meaningful number so long as the 
theory is conformal. Letting Ri be the radius (inversely proportional to the temperature) 
and i?3 be the radius, the infinite- volume limit i?i/i?3 — > is equivalent to the infinite- 
temperature limit in which Ri — > 0. In this sense we can compare our finite-volume results, 
at high-temperature, with known infinite- wo/ume, finite-temperature results. 

Introducing D7-branes orthogonal to the initial D3's in four directions introduces new open 
string degrees of freedom in the D3 worldvolume theory, from D3-D7 strings, and breaks half 
the super symmetry. We will denote the number of D7's as Nf. The so-called probe limit 
consists of keeping Nf fixed as oo. In this limit, the backreaction of the D7's on 

the geometry is negligible and hence the result in the near-horizon limit is Nf probe D7- 
branes embedded in AdS^ x S^. From the field theory perspective, knowing that the J\f = 4 
supersymmetry must be broken to J\f — 2 and from Chan-Paton factors of D3-D7 strings, 
one can argue that this corresponds to adding matter in the fundamental representation of 
the gauge group to the Af — A theory. Specifically, a number Nf of flavors in the form of 
M — 2 hypers have been added. 

Prom the field theory point of view, these hypers may be massive without breaking super- 
symmetry although of course this will break conformal invariance. This is manifest in the 
dual brane picture. In this setup, the D7 wraps all of AdS^^ as well as a trivial cycle inside 
the S^, namely an 5*^. A stable "slipping mode" then exists, which is simply a scalar field 
on the worldvolume of the D7 that depends only on the AdS^ radial coordinate [5] . The 
that the D7 wraps inside the may "slip off' the S^, as allowed by topology, i.e. contract 
to a point. If this occurs at a finite value of the radial coordinate, the hyper in the field 
theory can be shown to have a nonzero mass. The zero mass conformal case is then simply 
a D7 wrapping the same equatorial for all values of the radial coordinate. 

Indeed, via the AdS/CFT dictionary [5], this scalar slipping mode encodes both the mass 
of the hyper and the chiral condensate as a function of the mass, as we show below. The key 
observation is that the D7 worldvolume scalar has a mass (as a scalar in AdS^) of m? = —3 
(negative but above the Breitenlohner-Preedman bound [24]) and hence must be dual to a 
CFT operator of dimension 3 or 1, built from hyper fields. The only such unitary operator is 
a fermion bilinear. AdS/CFT then states that the source for the dual operator (the fermion 
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mass) and the value of the chiral condensate can be extracted from the asymptotic expansion 
of the scalar. These statements are made precise in section 3.3. 

Our goal is thus to map the phase diagram of this A/" = 2, non-conformal theory in the 
plane of mass versus temperature using this chiral condensate as our order parameter. In 
what sense is this chiral condensate an order parameter, however? That is, what are the 
symmetries of our theory and what symmetry breaking does this order parameter detect? 

Af — A SYM has an 5'C/(4) ~ 5*0 (6) R-symmetry. The hypers break the R-symmetry to 
SO{4:) X SO{2). This is easy to see on the gravity side: the D7 wraps an inside the S^, 
breaking the SO{Q) isometry to the S0{4) of the and an S0{2). This SO{2) ~ U{1) is 
a chiral symmetry in that the hyper fcrmions, which have opposite chirality, carry opposite 
charges under it. With A'^^ flavors, the theory has in total a global U{Nf) x S0{4)ji x 
S0{2)r ~ SU{Nf) X U{1)b X S0{4:)r x U{1)a. The U{1)a is our chiral symmetry in 
this case. At finite Nc, this U{1)a would be anomalous but as in large- A^'c QCD the axial 
anomaly is suppressed as Nc ^ oo. Without the Yukawa interactions required by J\f — 2 
supersymmetry the U{Nf) global symmetry would be enhanced to a U{Nf)L x U{Nf)R. 

In this paper, we consider only Nf — 1. This poses no problem on the gravity side: for 
this system, at zero mass, the factor of Nf would simply appear in front of the D7 action 
and would not affect the equations of motion^. These statements remain true at nonzero 
mass so long as all the flavor branes end at the same radius. A construction with different 
D7's ending at different radial values is possible, where the details of the spectrum would fix 
what remains of the U{Nf), but we will not consider such a situation here (though it would 
be a straightforward extension). 

3 The Probe Brane Action 
3.1 Fefferman- Graham Coordinates 

We will now write the action for the D7 slipping mode and show how to extract the mass 
and chiral condensate from solutions of the equation of motion (EOM). We also describe 
the boundary conditions in detail as these are somewhat different from those in [17, 18]. 
In these references and in much of the literature, a shooting method is employed, which 
is appropriate given only UV, or boundary, data. We impose boundary conditions in the 
IR, which is much simpler for numerics. In section 4 we will show that our IR boundary 
conditions reproduce the results of [17, 18] for the flat case. Boundary conditions similar to 
ours were used recently in [22] and also reproduced these results. 

In units where we set the curvature radius L to one the AdS^ metric in global coordinates 

^One subtelty that would however be present for Nf > 2 is that the configuration of lowest free energy 
might be found on the Higgs branch [19]. For Nf = 1 there is no Higgs branch since all hyper scalars have 
a quartic potential. 
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has the form 

ds" = GABdx^dx'' = TTT + fir)dT^ + r^dfll (3.2) 

fir) 

rf ^ _ f l+r'^ thermal AdS 

■^^"■^ ~ \ l + r^-^ AdS-Schwarzschild ^ ' 

dQp is the standard metric for ap-sphere. Here the time coordinate r has Euchdean signature 
and is periodic with period 27rT for temperature T. The parameter M is related to the black 
hole mass rribh as M = ^^rribh and to the temperature as M = nT. For AdS-Schwarzschild, 
the horizon is at r+, the largest solution to the equation /(r+) = 0. The flat case is the same 
with no 1 in /(r) and with Euclidean 3-space instead of the 5"^. 

In the above coordinates, the boundary is at r — > oo. This makes extracting asymptotics 
from numerics difficult. These coordinates are also not well-suited to holographic renormal- 
ization. We thus switch to a Fefferman- Graham [25] coordinate system, where the metric 
takes the form 

dz^ 1 

ds'^ = — + -^gij{z, X, M)dx'dx^ (3.4) 
z z 

gij{z,x,M)dx'dx^ ^ ^{1 - z\l + 4M'^)yF{z, M)-^dT^ + ^F{z, M)dfll (3.5) 

F(z,M) = 1 -2^2 + ^(1 + 4M^) (3.6) 

with radial coordinate z and x coordinates x*. The boundary is now at 2; = and the 
horizon is at 2;+ = (1 + 4M^)~^/^. Taking M — gives gij for thermal AdS with center at 
z — 1. Notice also that in our conventions our boundary has radius 1/2. 
The metric on the slice has an asymptotic expansion 



gi,{z,x,M) = g^ + g^^z' + gtjz' + Oiz') (3.7) 
1 

r 



4 = -dzag{l,S') (3.8) 



4 = -ldiag(-l,S') (3.9) 



gf = ]diag{l-12M\{l + 4M^)S^) (3.10) 
J 4 



where "yS"^" denotes the metric components. The metric is 

dnl = de^ + sin^ edijj^ + cos^ Odnj (3.11) 



10 



The D7 then wraps the whole of AdSr, and the inside the S^. The D7 shpping mode 
is the worldvolume scalar field 6 = 0{z), which is a function only of the radial coordinate to 
preserve the Lorentz invariance of the boundary theory. If ^(^o) = f for some finite Zq, then 
the has "slipped off" the S^. 

The D7 action is the Dirac-Born-Infeld (DBI) action, which in this background is simply 
the volume of the probe brane: 

Sd7 = Tnr J d^^^detP[G] (3.12) 

where P[G] is the puUback of the spacetime metric. With the 0{z) ansatz, the relevant 
part of the action is 

Sd7= [ dz^\F{z, M){1 - z\l + AM^)) cos^ 9{z)y/l + z^e'{z)^ (3.13) 
J 16 z^ 

Again thermal AdS is the same with M — 0. Our objective is to solve the 9{z) EOM. We 
show in section 3.3 how to extract the mass and chiral condensate from a solution and in 
section 3.4 we explain the boundary conditions we impose on solutions. 

3.2 Operator normalization 

While we have already identified the operator dual to the slipping mode 6 as the dimension 
3 fermion bilinear, there is still some freedom in the overall normalization. The correct nor- 
malization of the operator can be fixed by studying the 2-point function at zero temperature. 
Performing the integral over the in the 7-brane action gives a 5d action with prefactor 



Nf , 



where we used that in our L = 1 units (a') ^ = A = Oym^c = '^t^Qs^c- According to the 
AdS/CFT dictionary the dual operator O has a 2-pt function 

(O(.)O(0)> = J,H^_I(%-4 = ^J,. (3.15) 

^ ^ ' ^ 'I gl r(A - (i/2) At:^ ^ ^ 

In the field theory we are interested in the source term (the mass) and a vacuum expec- 
tation value for the operator O = AA where A is a canonically normalized fermion. So far 
we have only established that O = CO for some normalization constant C, but in order to 
interpret our results in the field theory wc need to know C. One way to fix C is to compare 
the two point functions of O and O. Since O is a protected operator (which follows from 
the non-renormalization of the flavor current 2-pt function established in [35] by supersym- 
metry), its 2-pt function can be evaluated in the free theory. For a canonically-normalized 
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fermion the position space propagator is simply (A(a;)A(0)) = times flavor and color 

delta functions. So in the free theory we then simply get from Wick contractions 

From which we can conclude that 

O^^O. (3.17) 

Prom this relation it becomes clear that all masses and condensates we find in this paper 
have to be multiplied by ^ to translate into field theory units. The phase transitions we 
study will happen at masses of order \/~\T. This is the natural scale since the mesons in 
the theory only have masses of order [26]. At zero temperature the mass we obtain this 
way also agrees with the energy of a straight string stretching from the horizon to the flavor 
brane. 

Note that the scaling of \ with Nf and Nc is expected simply from large Nc counting 
rules. The factor of A however is a little surprising at flrst sight. Besides its effect on the 
normalization of the operators it will for example lead to free energies proportional to A 
at large A, possibly indicating that in addition to the free value there is a 1-loop quantum 
correction but no further contributions from higher loops. 



3.3 Holographic Renormalization and Free Energy 

Given a solution of the EOM we need to extract the hyper mass and the chiral condensate. 
We use the method of holographic renormalization (which we call "holo-rg") to do this 
as well as compute the on-shell action. We will now briefly review the holo-rg procedure 
[27, 28, 29, 30]. 

The precise statement of AdS/CFT equates the lOD on-shell SUGRA action with the 
generating functional of the boundary field theory. The leading asymptotic value of a bulk 
field serves as a source for the corresponding field theory operator, i.e. derivatives of the 
SUGRA action w.r.t. leading asymptotic values will give boundary correlators. This will be 
made explicit for our 6{z) field below. 

Both the SUGRA action and the field theory generating functional are formally infinite: 
the SUGRA action suffers IR divergences while the field theory has UV divergences. Holo-rg 
proceeds first by regulating the SUGRA action. This means integrating not to the z = 
boundary but only to z = e. Local covariant counterterms are then added on the z = e slice, 
built from the metric induced on the slice, 7jj, curvature invariants thereof and fields on the 
slice, for instance our ^(e). The coefficients of these counterterms are fixed by requiring all 
divergences to cancel. Functional derivatives w.r.t. asymptotic values can then be taken. 
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followed by removal of the regulator, e — > 0, yielding renormalized boundary correlators. 
Covariance and all other symmetries can be maintained throughout this procedure. Holo-rg 
also gives us a way to compute a renormalized free energy of our M = 2 theory, since via 
AdS/CFT this is just equivalent to the on-shell SUGRA action. We find this much more 
elegant, in principle and in practice, than a background subtraction technique. 

Now to make all of this exphcit for our system. The holo-rg procedure for precisely this 
system, a probe D7 in AdS^xS^, was worked out in [31], so we quote those results. According 
to the AdS/CFT dictionary, the mass and condensate for the hyper fermions come from the 
asymptotic expansion 



9{z) = 9oz + e2Z^ + Q2zHogz + . . . (3.18) 

where the mass will be given by and the condensate will be fixed by a-nd d-^- In this 
case, the log coefficient 02 = -^RqOq [31] with Rq the Ricci scalar of g^. For us Rq — 24 so 
02 = 20q. The counterterms for this system are given in [31]^: 



Li = (3.19) 

^2 = -^^/lR, (3.20) 

Lz = -log(e)V7^(i?,,it:^^-^it:5 (3.21) 

^4 = (3.22) 

U = ^\ogd{e)^R,e{ef (3.23) 

We will denote the subtracted action as Ssub = Sdj + '^i Li and the renormalized action as 
Sren = liiiie->o Ssub- Notlcc that = and R^ = e'^Rg and to leading order ^(e) ~ OqS, 

so the last counterterm is finite. The first three counterterms arc needed to renormalize the 
volume of AdS^. The divergent piece of L3 has coefficient RijR^ — |i?o- 9ijj ^^i^ 

quantity vanishes and hence L3 is not needed. L4 and L5 are the counterterms for a free 
scalar in AdS^. The coefficient of is fixed by requiring a supersymmetric renormalization 
scheme [31]. Our L5 is actually slightly different from the one in [31], which has log(e) instead 
of log^(e). This is so a large- mass expansion correctly reproduces the flat case answer, as 
we now show. The condensate is given by [31] 

''This reference used a convention that AdS space had constant positive curvature. We use the opposite 
convention, hence our counterterms have —R^ relative to those in [31]. 
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= -2e2 + ie» + ^fl„e„iog(e2). (3.26) 

Let us focus on zero temperature. At large mass any finite volume effects should be 
neghgible and we should recover the flat space results, that is the leading piece in both vev 
and free energy vanish [31]. Writing the vev as a power series 

the coefficients oq, Oi, . . . can be solved for order by order in It is convenient to reinstate 
R3 in the bulk metric: 

,2 dz^ 1/1 2 1/1 2 z'^\ .^^o 

The equations of motion following from the DBI action 



can be solved order by order in powers of 



sin(6'(2;)) = argo{z) + -^argi{z) + .... (3.30) 

To leading order one reobtains the flat case solution argo{z) — cz and hence oq — since 
the vev in the flat case solution was zero. 
Expanding to order R^^ one obtains 

Iz-c^z^ + c^z^\og{zcf 
arg,{z) = ^-^^ . (3.31) 

First note that the prefactor of ^ compared to c in the 0-th order solution indeed nicely 
combines with the to form the dimensionless expansion parameter Also note that 

the solution contains a logc term, which nicely combines with the log 2; into logzc. It is 
straightforward to calculate {qq) for this analytic solution and we flnd with the counterterms 
above that oi = 0, so that in the large mass limit the vev vanishes^. It is this analytic answer 

''^Notc that 02 and all higher terms multiply a negative power of m . so that those contributions to the vev 
vanish at large mass even for non-zero coefficients. Indeed we have also worked out the 2nd order correction 
to 9{z) and found a non-zero 02. 
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for the first order correction that demands the particular form of L5 we use. (Note that this 
only works because the log terms combined in this particular way.) Any additional finite 
counterterm such as ^J^R.~^(P' would spoil the vanishing of the 02 term so all the counterterms 
are specified uniquely by the requirement that the vev vanishes in the large mass limit. It is 
easy to see that the analysis also carries through at nonzero temperature, where the metric 
only gets modified at order so that finite temperature does not affect Oq and 0\. 

With these counterterms implemented in the numerics the condensate at very large masses 
doesn't go exactly to zero and even dips negative. Our analytic solution makes it clear that 
this is a numerical artifact. It is easy to understand that the numerics becomes questionable 
at large mass, since in that case the position of the flavor brane pushes into the region in 
which we try to fit the numerical solution against the expected near boundary behavior in 
asymptotically AdS space. 

Even at small mass, however, we see the condensate becoming negative (at least for the 
curved case: see figures (6) and (8)). While at first sight this may seem disturbing, it is 
not forbidden. In non-relativistic quantum mechanics first-order perturbation theory always 
overestimates the free energy (ground state energy), i.e. the first-order correction is always 
positive. The second order correction is always negative and therefore the free energy is 
always a concave increasing function. A relativistic field theory, however, requires renormal- 
ization. Counterterms proportional to R^m? and will be present, and depending on the 
choice of their coefficients (choice of scheme) the corrections to the free energy may make 
the free energy decrease (and possibly convex). The condensate is the derivative of the free 
energy and hence may become negative. The scheme dependence of the field theory shows 
up holographically precisely in L5 and Lf. We believe we have fixed these counterterms in 
a physically well-motivated way and that negative condensates are simply a consequence of 
this choice of scheme and no cause for alarm. 

The free energy of the field theory is given by the on-shell value of the renormalized 
action. The order A'^^ contribution to the free energy comes from the lOd SUGRA action 
and is untouched by the introduction of the probe brane. Naively there are two sources of 
order NfNc corrections, the renormalized value of the on-shell DBI action as well as the order 
^ correction to the order on-shell bulk SUGRA action due to the backreaction. The 
latter however is proportional to the variation of the supergravity action with respect to the 
unperturbed background metric and hence vanishes. The backreaction only enters at order 
N'^ and the only order NfNc contribution to the free energy comes from the renormalized 
on-shell DBI action. 

3.4 Boundary Conditions 

We will discuss first the M = thermal AdS case, with no black hole horizon. We consider 
two kinds of solutions. The first kind are D7's ending at some finite zq, i.e. ^(^o) — 7r/2. 
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This is our first boundary condition, and we want a second condition in the IR. We find one 
by requiring that the D7 ends smoothly, with no conical deficit. The pulled-back D7 metric 
contains terms 

d^m ^ (4 + G\^?)dz^ + cos^ e{z)dnl (3.32) 
z 

Let a — cos9{z) so that da — —sm9{z)9'{z) and hence 

dsl,D^-p-^da' + a'dnl (3.33) 
0'[zy sm 

No conical deficit means the coefficient of da'^ must be one when z = Zq- This gives 
Zq'^9'{zq)^'^ = and hence 9'{zo) = oo. We implement this numerically with a number on 
the order of 10^. 

The second kind of solution is for D7's extending to the center of AdS at 2:0 = 1. In this 
case, we choose a boundary condition 9{z — 1) = /3 for some angle 0</3<7r/2. P — 
is the massless case and dialing P will interpolate smoothly to a D7 ending precisely at the 
center. The second boundary condition is now the Neumann condition 9'{z — 1) = 0. 

In AdS-Schwarzschild we also have two kinds of solutions, D7's ending outside the horizon 
or ending on the horizon. For D7's ending outside the horizon we use the same boundary 
conditions as those in thermal AdS. For D7's ending on the horizon we use boundary con- 
ditions 9{z+) = (3 and 9'{z+) = 0. The fiat case, where the horizon is present for any finite 
temperature, also has these two types of solutions with these boundary conditions. 

In summary, our procedure is this: solve the EOM for 9{z) with these boundary conditions. 
With this solution, we can immediately find the mass and condensate for that value of zq or 
p. We also plug this solution into Sren to find the free energy. We can then go to the next 
value of Zq or f3, which amounts to dialing the mass. In this fashion, from the behavior of the 
condensate and the free energy, we can identify the chiral phase transition and determine its 
order. 

In much of the literature (for example [17, 18, 32]), IR boundary information such as this 
is not used. Instead, a shooting method is used. To justify and cross-check our method, in 
the next section we will reproduce some known results for the flat case using our boundary 
conditions and holo-rg. We will also compare our curved case results in the high-temperature 
hmit to these results as a useful check. 



4 The Flat Case 

The thermodynamics of this J\f — 2 theory has been studied in [17, 18]. We will compute 
the chiral condensate as a function of the mass, the free energy as a function of the mass 
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vev vev 




Figure 3: (a.) Chiral Condensate as a function of mass (M = 1); gray and black curves 
indicate solutions ending inside and outside the horizon respectively (b.) Close-up view of 
(a.) 

and then draw the phase diagram. As suggested in [17, 18], a first-order phase transition 
does occur at a critical value of the mass. The critical value is merit ~ 1-30 (for M — 1)^. 
The form of the free energy shows this clearly. In this case, our AdSr, metric is 

^ ^ + li(l - z'M^Fiz, MY'dr' + l^F(z, M)dx^ (4.34) 
z-^ z^ I z^ I 

F(z, M) = 1 + z^M^ (4.35) 

and the horizon is at z+ — 1/M. An important fact we will need later is that the last term 
in the metric has a 1/2 whereas the curved case metric had a 1/4, hence z in the flat case 
is ~ \/2z in the curved case and the mass we extract from eq. 3.18 in the flat case will thus 
be 1/ \/2 times the mass we compute in the curved case. This is important when comparing 
the phase diagrams in the two cases, in particular their slope, as we do in section 5. 
The D7 action is now 

Sd7 = / dz~{l - z^M^) cog^9{z)^l + z''e'{zY. (4.36) 

J 4:Z^ 

The counterterms are given above but with L2 — — — because — 0. 
Our results for the condensate are summarized in flgure (3) and for the free energy in (4), 
clearly reproducing the first order phase transition found in earlier work. The condensate as 

^Our merit is larger than the merit ~ 0.92 in [17, 18, 22] by a factor of \/2. This comes from our choice 
of coordinates. The coordinate r in [22] is related to onr z by z = -^i. We verified that this change of 
coordinates reproduces fig. 2(b) of [22]. While a change of coordinates clearly cannot change a physical 
quantity like the mass, different coordinate systems suggest different natural defining functions to extract 
the boimdary metric, in our case but for [22]. The corresponding boundary metrics differ by a factor of 
V2 for the radius of the 5*"^ and hence the mass measured in units of the inverse radius of the sphere differs 
by a factor of as well. 
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(a) (b) 
Figure 4: (a.) Free Energy as a function of mass (M = 1). (b.) Close-up view of (a.) 
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Figure 5: Phase Diagram of the N — 2 theory in fiat space. T — Mjii. 

a function of mass is multivalued; the plot of the free energy shows that at a critical value 
of the mass we discontinuously jump from a D7 ending some finite distance away from the 
horizon to a D7 touching the horizon. In terms of the induced metric on the world volume 
of the D7 brane this transition corresponds to a topology-changing phase transition. 

The phase diagram for this theory appears in figure (5). As explained in the introduction, 
this is a straight line. A linear fit to this data gives a slope of ~ 4.1 and a y-intercept of 
zero. This is expected since in a conformal theory two scales are needed for a transition and 
either T or m going to zero eliminates one scale. 



5 The Curved Case 

We start with thermal AdS, for tempratures below T^p. We find a surprise. In terms of the 
topology of the D7 brane there are still two classes of solutions: the D7 brane can either end 
at a finite value of the radial coordinate and hence have a finite size S"^ inside the AdSs, or 
it can reach all the way to the center of AdSs with this time the S'^ inside the S'^ remaining 
finite. The two branches meet at the point where both 5'^'s shrink at the center of AdS. In 
the field theory we expect the two configurations to correspond to a competition between 
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Figure 6: (a.) Thermal AdS chiral condensate as a function of mass. Gray curves are branes 
ending at the center of AdS, black is ending at finite z < 1 (b.) Close-up of (a.) 
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(a) (b) 
Figure 7: (a.) Thermal AdS free energy as a function of mass, (b.) Close-up of (a.) 



the explicit mass terms we added for the hypers with the curvature induced scalar masses. 
Figure (6) shows our results for the vev and figure (7) shows the free energy. Neither the 
free energy nor the condensate is multiple-valued in the mass. In other words, despite the 
topology-changing transition in the bulk we find no first-order transition in the field theory. 

Note that our free energy has a non- vanishing positive value of 0.09 ~ 3/32 at zero 
temperature. This can be interpreted as the XNfNc contribution to the Casimir energy of 
the system on the sphere. In the massless case the solution can even be obtained analytically, 
9{z) — 0, and we get an analytic expression for the DBI contribution to the Casimir energy, 

where V denotes the volume of the 3-sphere. The leading order (in Nc) contribution to the 
Casimir energy comes purely from the supergravity action evaluated on background AdS5 x 
geometry itself and has been obtained in [29] , 
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(5.38) 



Their result matches perfectly the free field theory result, where one basically just counts 
the number of fields and hence gets a contribution of order from the adjoint fields with 
no powers of A. This agreement is guaranteed by supersymmetry. In [33, 34] it was shown 
that the Casimir energy can be obtained from the conformal anomaly by exploiting the fact 
that the sphere is conformal to flat space and the latter is assigned zero Casimir energy. The 
conformal anomaly however is completely protected (that is, independent of the coupling 
constant) in a theory with at least M = 2 supersymmetry^ and hence so is the Casimir 
energy on the sphere (or any other space conformal to fiat space). 

The order NfNc correction we found to this leading order Casimir energy is proportional 
to A and hence can not possibly match the free field value, but rather seems to come from 
a 1-loop contribution with no higher-loop corrections. Since we have J\f — 2 supersymmetry 
the conformal anomaly should still be protected^. However, this time there is an additional 
contribution to (T^) from the scale anomaly due to the non-vanishing beta function. In an 
M = 2 theory the beta-function is 1-loop exact and hence sees only the order A correction 
even at strong coupling. Our result suggests that the Casimir energy might in a similar 
fashion only get a contribution from 1-loop. It would be interesting to see if this can be 
verified from a weak coupling calculation. 

Finally we come to the AdS-Schwarzschild case ior T > Thp- We find that for T > 4*Tffp 
the free energy and condensate already have nearly the same form as the fiat case in figures 
(3) and (4). As T decreases, these continuously deform to the shapes shown in figures (8) and 
(9) where T = 1.75 * Thp-^ In close-up we flnd that the condensate is again double-valued 
and will jump discontinuously. The free energy again has a discontinuous flrst derivative just 
as in the flat case. This remains true all the way down to T^p so no critical point analogous 
to the point C in the large- A^c QCD phase diagram appears for this theory. 

Putting both the high temperature and low temperature regions together we finally arrive 
at the phase diagram for our theory, as displayed in figure (2). No first-order transition 
occurs below the T = Thp vertical dcconflncment transition line. A linear fit indicates that 
the chiral transition line has slope ~ 5.8. This agrees with the earlier result of 4.1 from the 

^With Af = 1, supersymmetry still relates the conformal anomaly to a protected R-current anomaly, but 
different U(l)s can play the role of the superconformal R-current at different values of the coupling. 

^Indeed it was shown in [35] that in a closely related system with an 07 in addition to 4 probe D7s 
the supcrgravity anomaly perfectly reproduces the field theory anomaly. It is easy to sec that without the 
orientifold their analysis for a single D7 reproduces the right value corresponding to a fundamental hyper. 
In their case the field theory was exactly conformal (not just to leading order in ^) and in the bulk the 
order XNfN^ contribution to the Casimir energy canceled between the orientifold and D7s 

^At small mass the condensate begins to dip to negative values around 3.7 * Thp- It may also dip 
negative at higher temperatures, but our numerics could not achieve sufficiently high resolution to see this 
for T > 3.7* Thp. 
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(a) 

Figure 8: (a.) AdS-Schwarzschild chiral condensate as a function of mass. Gray curves are 
branes falling into the horizon, black is ending outside. T — 1.75 *Thp (b.) Close-up of (a.) 
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flat case after taking into account the factor of \/2 difference between the z coordinate of 
our flat-case metric and that of our AdS-Schwarzschild metric: 5.8 ~ \/2 * 4.1. The chiral 
transition hne intercepts the deconflnement hne at m ~ 2.6. 

6 Conclusion 

We have presented a numerically-generated phase diagram for large- A^c A/" = 4 super- Yang- 
Mills on a 3-sphere coupled to N j massive fundamental hypers, including the first-order 
chiral phase transition line. Below the Hawking-Page temperature we find that, despite 
a topology-changing transition in the bulk, a first-order transition does not occur in the 
boundary theory. 

Isospin chemical potential could be included in this analysis by introducing multiple D7's 
ending at different values of the AdS radius, allowing for 7-7 strings and hence heavy-light 
mesons. Baryon number chemical potential could be included by turning on the gauge po- 
tential on the D7 worldvolume, although it is unclear what boundary conditions to impose on 
this field in the IR. Supposing that question could be answered, this would be one advantage 
AdS/CFT would have over lattice gauge theory which currently cannot easily include baryon 
number chemical potential. Other interesting extensions include studying the corresponding 
phase diagram at weak coupling and studying the effects of the D7 backreaction. 
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